Plaquette-RVB state in the Frustrated Honeycomb Antiferromagnet 
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We study the proposed plaquette-RVB (pRVB) state in the honeycomb lattice Ji — J2 model with 
frustration arising from next-nearest neighbour interactions. Starting with the limit of decoupled 
hexagons, we develop a plaquette operator approach to describe the pRVB state and its low energy 
excitations. Our calculation clarifies that the putative pRVB state necessarily has f-wave symmetry 
- the plaquette wavefunction is an antisymmetric combination of the Kekule structures. We estimate 
the plaquette ordering amplitude, ground state energy and spin gap as a function of J2/J1. The 
pRVB state is most stable around J2/J1 ~ 0.25. We identify the wavevectors of the lowest triplet 
excitations, which can be verified using exact diagonalization or DMRG studies. When J2 is reduced, 
we can have either a deconfined Quantum Phase Transition (QPT) or a first-order transition into 
a Neel state. When J2 is increased, we surmise that the system will undergo a first order phase 
transition into a state which breaks lattice rotational symmetry. 
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I. INTRODUCTION 

Magnetic systems with frustrating interactions serve 
as an excellent route to generating novel quantum states 
with various degrees of entanglement. They often give 
rise to many body ground states in which many con- 
stituents become entangled, from which new collective 
degrees of freedom emerge. Examples include Valence 
Bond Solid (VBS) phases [ll, decoupled chainsQ and spin 
liquids. At each level of entanglement, a new effective 
theory is needed to describe the low energy properties of 
the system. In this light, we study a plaquette-ordered 
phase proposed as a ground state for the honeycomb lat- 
tice spin-1/2 Ji — Ji model, with frustration arising from 
next-nearest neighbour couplings. This plaquette-RVB 
(pRVB) phase comprises of entangled hexagons, arranged 
m a pattern. 

Recently, the J1-J2 model on the honeycomb lattice 
has generated tremendous interest as an effective model 
for the intermediate-U physics of the Hubbard model on 
the honeycomb lattice [3(. It may also be relevant to the 
material Bi3Mn40i2(N03) which shows no long range 
order down to the lowest temperatures [3-0] ■ Being the 
simplest model of frustration on the honeycomb lattice, 
the J\ — J2 model has been extensively studied. In the 
semiclassical limit of large-S, this model is very well un- 
derstood: for small J 2 , the ground state has Neel order. 
At a critical value J 2 = 1/6 (we set Ji = 1), there is a Lif- 
shitz transition into a spiral state with extensive ground 
state degeneracy. Order-by-disorder effects mediated by 
quantum/thermal fluctuations select three spiral states 
which are related by lattice rotations. The ground state 
thus breaks lattice rotational symmetry. By contrast, 
the extreme quantum limit of S — 1 /2 is not well under- 
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stood. In spite of many recent studies, the phase diagram 
has not yet been conclusively established. In particular, 
many studies point to an interesting intermediate phase 
for 0.2 <J 2 <0.4 which interpolates between Neel order 
and a state which breaks lattice rotational symmetry. 

While all studies in the literature agree upon the ex- 
istence of an intermediate phase, its precise nature has 
not been established. Some proposals indicate that the 
intermediate state may be a i? 2 spin 

liquid|M3- Oth- 
ers give support to apRVB state which breaks transla- 
tional symmetry [llTfl3j . A recent variational calculation 
argues that the intermediate state does not have plaque- 
tte order[14J. In the light of such conflicting studies, we 
study the proposed pRVB state and suggest more precise 
tests to confirm its existence. It is worthwhile to devise 
such tests as most of the evidence for a pRVB ground 
state comes from calculations restricted to small system 
sizes. In addition, our study sheds light on the properties 
of this putative state, its low energy excitations and the 
nature of phase transitions into and out of this phase. 

In this paper, we build an effective theory for the can- 
didate pRVB state using a plaquette- operator approach. 
We estimate its ground state energy and spin gap and de- 
termine the nature of low-energy excitations. Our goals 
are threefold: (i) to understand the nature of the pRVB 
state, (ii) to clarify the nature of phase transitions into 
and out of this state, and (ii) to suggest tests which can 
be carried out using numerical techniques such as the 
Density Matrix Renormalization Group (DMRG). Such 
tests might conclusively reveal whether the pRVB state 
indeed arises in the J1-J2 model. 

This paper is organized as follows. In Sec. [Til we begin 
with the S=l/2 J\ — J 2 model on a single hexagon. We 
classify all its eigenstates as a function of J2/J1. Sec. 
IIIII introduces the plaquette operator representation and 
lays out the formalism sequentially - by first consider- 
ing two coupled hexagons and then tiling hexagons to 
form a honeycomb lattice. We argue that the pRVB 
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FIG. 1: The proposed pRVB state on the honeycomb lattice. 
The shaded hexagons are in an anti-symmetric combination 
of the two singlet covers 



state has f-wave character, rather than s-wave. We ar- 
rive at a simplified effective Hamiltonian that only in- 
volves triplet excitations. Sec. [TV] gives the results of the 
pRVB mean-field theory with spin-gap, ground state en- 
ergy and pRVB amplitude. Finally, Sec |V] discusses the 
nature of phase transitions into and out of this phase and 
implications for the Ji — J2 model. 

II. THE SINGLE HEXAGON PROBLEM 

On a single hexagon, the spin-1/2 J1-J2 model is easily 
diagonalized, with a Hilbert space of 64 (= 2 6 ) states. 
The Hamiltonian is given by (we set Ji=l) 

H s ingle~hex = ^ Sj • Sj + Jg S Si ■ Sj. (1) 

We take the exchange couplings Ji and J2 to be positive. 
The Ji and J2 bonds are shown in FigHJ 

As the Hamiltonian is invariant under global spin rota- 
tions, its eigenstates have well-defined magnetic quantum 
numbers. Diagonalizing the Hamiltonian numerically, we 
obtain the spectrum consisting of 5 singlets (St o t=0), 
27 triplets (S tot =l), 25 quintets (S tot =2) and 7 heptets 
(S tot =3). In the next section, we will argue that only one 
singlet and the 27 triplet states are relevant for under- 
standing the pRVB phase. 

The hamiltonian is invariant under two spatial opera- 
tions: (i) rotation by 60°, and (ii) mirror inversion about 
a line (see FigJ5|). As the two operations do not com- 
mute with each other, the eigenstates of the hamiltonian 
cannot be chosen to be simultaneously the eigenstates of 
both. For later convenience, we choose the eigenstates to 
be simultaneously eigenstates of rotation. Defining R to 
be the rotation operator, we note that R is the identity 



operation. Therefore, the eigenvalues of R are the sixth 
roots of unity. 

We denote the states obtained by diagonalization as 
I u) , with u = 1 , . . . , 64 with the energy eigenvalues 
e u . Each state is assigned three quantum numbers 
(l u , m u , f u ) where l u , m u are the usual magnetic quan- 
tum numbers and r u is the eigenvalue under rotation. 
For example, {l u ,m u ,r u ) — (1,0,-1) indicates that the 
state I it) is a triplet with S z — and it changes sign 
under rotation by 60°. Note that these three quantum 
numbers do not uniquely identify a state. 

As J 2 is varied, some of the the low energy states of 
the single hexagon problem are: 

(i) '+' singlet: this state is the symmetric combination 
of the two Kekule structures proposed for benzene. With 
quantum numbers (0, 0, 1), it is invariant under rotation 
by 60°. 

(ii) '-' singlet: this is the ground state for J2 < 0.5. 
It is predominantly (~ 98.5% - the precise fraction de- 
pends weakly on J2) an antisymmetric combination of 
the Kekule structures (see FigllJ. With quantum num- 
bers (0,0, —1), it is even under inversion and odd under 
rotation. 

(iii) t\ triplets: these are the triplet states with the 
quantum numbers (1, m, 1), with m = — 1, 0, 1. They are 
invariant under rotation. 

Fig [5] shows the low-energy spectrum of the single 
hexagon problem as a function of Ji- We only indicate 
the ground state and the first two excited states. The 
'-' singlet is the ground state for J2 < 0.5. There is a 
level crossing at J2 = 0.5 beyond which the '+' singlet 
becomes the ground state. 

On the full honeycomb lattice, previous studies have 
suggested a plaquette-ordered ground state for 0.2 < 
Ji J$ 0.4 |14] - [l3j . In this parameter regime, the ground 
state of the single hexagon problem is the '-' singlet. 
We argue that the suggested plaquette-ordered state is 
composed of a V3 X \fi ordering of hexagons in the '-' 
state as shown in Fig[TJ The choice of the '-' state corre- 
sponds to f-wave pRVB order - the plaquette wavefunc- 
tion changes sign upon rotation by 60°. In the next sec- 
tion, we describe the plaquette-operator method which 
provides a consistent scheme to include quantum fluctu- 
ations around this state. 



III. PLAQUETTE OPERATOR THEORY 

We obtain an effective theory of the pRVB state in 
analogy with the bond operator formalism of Sachdev 
and BhattJUj]. Developed in the context of dimer states 
in spin-1/2 systems, the bond operator method starts 
from the limit of isolated bonds - the ground state is then 
composed of disconnected dimers (singlets). A bosonic 
representation is introduced for the states of each bond 
in the singlet-triplet basis, with a constraint of unit boson 
occupancy to preserve the Hilbert space. The dimer state 
can now be elegantly described as a Bose-condensate of 
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FIG. 2: Top: Exchange couplings in the single hexagon 
Hamiltonian. The following operations commute with the 
Hamiltonian: (bottom left) rotation by 60 °, and (bottom 
right) inversion about an axis passing through opposite sites. 



We call these 'plaquette-bosons'. To preserve the Hilbert 
space, we impose a constraint of unit boson occupany. 
On a single hexagon, we may write 



bt\Q). 



(2) 



Here, the operator b^ u creates a u-boson when it acts upon 
the vacuum. The vacuum itself is an unphysical state, as 
it does not respect the unit occupancy constraint. Oper- 
ators can also be translated into the new representation, 
e.g., \v)(u\ = b\b u . 

Having solved the single hexagon problem exactly, we 
tile hexagons to form a two-dimensional honeycomb lat- 
tice, as shown in Fig. [TJ The shaded hexagons serve 
as sites upon which plaquette-bosons reside. As a start- 
ing point, we consider the limit of isolated plaquettes. 
Turning off all inter-plaquette couplings, we obtain a 
quadratic theory of bosons given by 



H 



intra— plaq. 



(3) 
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FIG. 3: Energy spectrum of the single hexagon problem - for 
brevity, we only show the lowest three states. The ground 
state is the '-' singlet for J2 < 0.5. For Ji > 0.5, the '+' 
singlet becomes the ground state. For J2 < 0.25, the first 
excited states are the t\ triplets. For J2 > 0.25, the ti go 
higher in energy. 



singlet bosons. Moving away from the limit of isolated 
bonds, inter-bond coupling can be naturally interpreted 
as four-particle interactions between bosons. Assuming 
that dimer order is strong, one arrives at a mean- field pre- 
scription to treat the interaction terms which gives rise 
to a quadratic theory of 'triplons'. Corrections beyond 
mean- field theory, viz., interactions between triplons, can 
be ignored as the triplons are present in dilute concentra- 
tions. This theory has been very succesful in obtaining 
the spectrum of triplon excitations above the dimer state. 
The dimer-Neel phase transition is easily understood as 
the Bose-condensation of triplons with the bond-operator 
mean-field theory providing excellent quantitative esti- 
mates for the transition point. 

This approach has been adapted to study plaque- 
tte order, albeit in a simple system with four-site 
plaquettes [T|| . Here, we develop a plaquette operator 
approach to describe six-site plaquette order in the hon- 
eycomb pRVB state. 

Having obtained all eigenstates of the single hexagon 
problem in the previous section, we now introduce a 
bosonic representation for these states. We associate a 
boson with each of the 64 states in the Hilbert space. 



where the index i sums over all shaded plaquettes in FigQ] 
- note that the shaded hexagons form a triangular lattice. 
The index u sums over the 64 eigenstates of the single 
hexagon problem. The operator 6 J creates a u boson at 
the site (plaquette) i. 

The pRVB state can now be identified as a Bose- 
condensate of the '-' bosons; i.e., we may set 



= u1 - 



P- 



(4) 



where p is the Bose-condensation amplitude. Physically, 
the quantity p 2 gives the probability of finding a shaded 
plaquette in the '-' state. We take p to be real as its phase 
is not relevant in the present context. In the isolated 
plaquette limit, p 2 is unity as every hexagon is entirely 
in the '-' state. 

To preserve the Hilbert space, the total occupancy of 
bosons per plaquette should be unity. In the spirit of 
mean-field theory, we enforce this constraint on average 
using a chemical potential /i. 



H -> H - fi 



Np 2 



E 



h \.uKu 



N 



(•5) 



,64 



where N is the total number of shaded plaquettes in 
FigQ] Eventually, we will set fi by demanding 



d(H) 



= 0, 



(6) 



which will tune the average boson occupancy to unity. 

The next step towards obtaining a consistent theory is 
to introduce inter-plaquette terms. We deduce the form 
of these terms by considering two coupled hexagons in the 
following section. Before giving details of our plaquette 
operator theory, we illustrate the approach by an analogy 
to spin wave theory. Consider a spin-S system that has 
Neel order at zero temperature. To zeroth order, this is 
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FIG. 4: Coupling between two 'nearest-neighbour' hexagons, 
A and B. The dotted lines denoted Heisenberg-like exchange 
terms which couple A and B. 



described as a classical system of anti-aligned spins with 
each spin of length S. This is analogous to our isolated 
plaquette limit with p = 1 which has perfect plaquette 
ordering. Standard linear spin wave theory takes into 
account quantum fluctuations about the Neel state and 
gives the spectrum of low-lying excitations. Spin waves 
reduce the length of the ordered moment and give rise 
to zero-point energy corrections. Analogously, our pla- 
quette operator theory captures the low-lying excitations 
about the pRVB state. These corrections will reduce the 
value of p and modify the ground state energy in a self- 
consistent manner. 



A. Coupling two hexagons 

FigQ] shows two nearest-neighbour hexagons, A and 
B, and the five bonds which couple them. The five 
Heisenberg-like bonds act as a perturbation, which we de- 
note as Hab , coupling the single-hexagon states of A and 
B. In the absence of this perturbation, the total eigen- 
state is a direct product of single hexagon states centred 
on A and B - \A U B V ), signifying that the A hexagon is 
in the ip u state and B is in the ip v state. 

The matrix element of Hab between two direct prod- 
uct states, \A U B V ) and \A p B q ), is given by 

(A u B v \H AB \A p B q ) = 
]T MAuB v \S lb , A \A p B q ). (7) 

b a—x,y,z 

The index b sums over the inter-hexagon bonds which 
connect site ib of hexagon A and site jb of hexagon B (see 
Fig. 0}. Each bond has an exchange coupling strength Jb 
(given by J\ for one bond and Ji for four other bonds). 
This quantity reduces to a product of matrix elements 



of single-hexagon operators, 
(A u B v \H AB \A p B q ) = J2 M A u\ s h,A,a\ A p) 

b a—x,y,z 

x(B v \S jb , B , a \B q ). (8) 

Using EqlH we arrive at the following lemma. Note 
that our motivation is to study the pRVB state in which 
each hexagon is predominantly in the '-' state. 

Lemma: If hexagon A is initially in a singlet state, the 
coupling hamiltonian Hab will take it to a triplet state. 
In particular, if A is initially in the the '-' state, then the 
coupling matrix element in Eq(5]can be non-zero only if 
A u - the final state - is a triplet. This is a consequence of 
the rotation properties of the spin operator Si b .A, a which 
acts as a triplet under rotation. 

We make the following assertions using this lemma. 

(i) (A-B-\H AB \A-B-) = 0. Thus, the strict plaque- 
tte RVB state which has both hexagons in the pure '-' 
state is not preserved under the coupling Hamiltonian. In 
the bosonic representation, the inter-plaquette coupling 
does not lead to terms of the form b A _b B _bA,-bs.- 
(equivalently, p A ). In the condensed state, after taking 
inter-plaquette terms into account, there is no term pro- 
portional to p 4 in the Hamiltonian. 

(ii) Similarly, (A-B-\H AB \A-B U ) = 0. There are no 
terms proportional to p 3 in the Hamiltonian. 

(hi) (A U B V \H ab\A-B _) can be non-zero only when u 
and v are triplet states. Furthermore, from conservation 
of z-component of spin, we should have m u + m v = 0. 
These matrix elements lead to 0(p 2 ) terms which involve 
creation/annihilation of pairs of triplets, i.e., we have 

terms of the form p 2 {b\ u b^ v + 

(iv) (A-B v \Hab\A u B-) can be non-zero only when 
u and v are triplet states. This gives rise to 0(p 2 ) 
terms which involve hopping of triplet operators, i.e., 

p 2 {&],Au}- 

(v) We can have non-zero matrix elements of the form 
(A_B u \HAB\A v B q ) - leading to 0(p 1 ) terms. However, 
these processes are highly unlikely and can be ignored. 
This is justified when the pRVB order is strong and each 
plaquette to be predominantly in the '-' state. While the 
'-' boson is condensed, all other bosons will be present in 
very dilute concentrations. Thus, interactions between 
such bosons are unlikely. 

(vi) We also have terms independent of p which do not 
involve the '-' state, but these can again be ignored using 
the above argument. 

The above statements drastically simplify the effec- 
tive coupling between adjacent plaquettes A and B. Our 
scheme can be justified in two different ways. In the first 
approach, we expand the hamiltonian in powers of p. Ar- 
guing from (i) through (iv) above, we see that the lead- 
ing terms are 0(p 2 ). These terms are quadratic in triplet 
operators and involve hopping-like and pairing-like terms 
discussed in (iii) and (iv) respectively. We retain these 
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leading terms alone and ignore 0(p),0(l) corrections as 
discussed in (v), (vi) above. 

A second approach is to start with the assumption of 
strong pRVB order. Every hexagon is predominantly oc- 
cupied by the '-' boson. All other bosons are present in 
very dilute quantities. Thus, beginning with (v) and (vi) 
above, we neglect unlikely interactions between these di- 
lute bosons. Using (i)-(iii), we find that the most likely 
processes are 0{p 2 ) terms which are quadratic in triplet 
operators. 

The above arguments apply to any pair of nearest- 
neighbour plaquettes in the honeycomb lattice of Fig|T] 
They lead to an effective hamiltonian for the pRVB state 
which only involves the '-' singlet and all triplet bosons. 
The quintets, heptets and other singlets simply drop out 
of the problem. Interaction corrections arising from 0(p) 
terms that we have ignored can populate these states. 
However, we expect these corrections to be small. In the 
rest of this paper, we ignore the other bosons as they 
play no role. 



B. Effective Hamiltonian for the pRVB state 

We have examined the nature of the coupling between 
nearest neighbour plaquettes, and deduced that they lead 
to hopping-likc and pairing-like terms in triplon opera- 
tors. There arc no couplings between further neighbours 
as the microsopic spin Hamiltonian has only short range 
interactions. A cursory inspection of Fig. [T] shows that 
Ji and J2 interactions do not couple further neighbour 
plaquettes. 

The effective theory of the pRVB state can be written 

as 

# = E (eu-/i)6.U^u + A r p 2 (e_- M ) + 7V M 

i u— triplets 

+p 2 J2 E [ b Uhv + (<e) b{ u b] tV + h. c ] .(9) 

(ij) u,v— triplets 

The first line includes the intra-plaquette terms, the in- 
dex i summing over all shaded plaquettes in Fig HI The 
chemical potential [i can be tuned to enforce single bo- 
son occupancy (on average). N denotes the number of 
shaded plaquettes in the system. 

The second line involves inter-plaquette terms - the in- 
dices u and v sum over all 27 triplet states in the single- 
hexagon spectrum. As argued at the end of the previ- 
ous section, the quintets, heptets and singlets other than 
the '-' state drop out of the problem. The coefficients 
(huv) an d (d l J v ) are the amplitudes for triplet hopping 
and triplet pairing processes. The evaluation of these 
amplitudes is discussed in the Appendix. 

We rewrite the Hamiltonian of Eq(3] by going to mo- 
mentum spacefill- As the plaquettes form a triangu- 
lar lattice, we obtain a hexagonal Brillouin zone. Fig. 
[5] shows the plaquette-Brillouin zone in relation to the 
hexagonal Brillouin zone of the underlying honeycomb 



lattice. The Hamiltonian takes the form 
T 



( b U \ 

"27, k 

h, -k 
V 627 ,-k J 



fle lMl fc(k) H d (k) 



027, k 



(10) 



The matrix H e ^.h involves triplet energies of the single 
hexagon problem, the chemical potential /1 and hopping 
terms h l J v . The off-diagonal blocks Hd involves pairing 
amplitudes d^ v . The primed sum indicates that if k is 
included, — k is to excluded. 

The eigenvalue spectrum of this hamiltonian can be 
obtained by a bosonic Bogoliubov transformation. These 
give the energies of 'triplon' excitations above the pRVB 
state. The ground state energy acquires zero-point en- 
ergy contributions from these triplon modes. The hamil- 
tonian has two parameters p and /1 which are determined 
self-consistently, by demanding 



d(H) 

d(H) 
dp 



0, 



= 0. 



(11) 

(12) 



EafTTI tunes the average boson density. We treat p as 
a variational parameter - Eq. Q2] minimizes the ground 
state with respect to p. 

For any value of J2 < 0.5, the hamiltonian in Eq|H]can 
be derived and diagonalized. (For J2 > 0.5, our plaquette 
operator theory is not justified as the '-' state is not the 
ground state of the single hexagon problem.) EqsQTJ[T2] 
can then be used to obtain the self-consistent solution. In 
the self-consistent solution, the eigenvalue of the lowest 
triplon mode gives the spin gap above the pRVB state. 



IV. RESULTS 

The plaquette operator mean-field theory is a novel 
method to study plaquette-ordered phases. Our plaque- 
tte operator calculation is valid in the limit of strong 
plaquette order, i.e., when p 2 is close to unity - as p 2 
is the probability of finding a plaquette in the '-' state. 
As a consistency check, we plot the self-consistently ob- 
tained value of p 2 as a function of J2/J1 in FigJSJ We find 
that pRVB order is strongest around J 2 ~ 0.25, where 
p 2 ~ 96%. As seen in the figure, we find p 2 > 90% 
for 0.15 < J2 < 0.43. This puts our plaquette operator 
approach on a solid basis. 

In FigJHl we plot the ground state energy obtained from 
plaquette operator theory as a function of J2/J1. This 
is obtained as the zero temperature expectation value of 
the Hamiltonian in Eq. [TUJ To see if the pRVB state is 
a plausible ground state, we compare its energy with the 
following estimates: 



6 



0.96 p- , , , , , , , , 

0.95 - ^ \ 

0.94 / \ 

0.93 / \ 

0.92 \ 

0.91 - \ 

0.90 \ 

0.89 \ 

0.88 I—. > ' ' ' ' ' — 

0.15 0.20 0.25 0.30 0.35 0.40 0.45 

J 2 

FIG. 5: The probability of finding a shaded plaquette in the 
'-' state, p 2 as a function of J2/J1. pRVB order is strongest 
around J2 ~ 0.25. 
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FIG. 7: The spin gap obtained from plaquette operator theory 
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FIG. 6: The energy of the pRVB state obtained from plaque- 
tte operator theory. We compare energy per unit cell with 
other candidate states (see text). 

(i) Neel/Spiral state from large-S spin wave theory: 
in the semi-classical large-S limit, the ground state has 
Neel order for J 2 < J1/6 and spiral order for J 2 > J1/6. 
Linear spin wave theory gives the energy of low lying 
excitations [18]. Fig|6] plots the energy of the large-S 
ground state including zero-point contributions from spin 
wave excitations. To obtain the energy for the spin-1/2 
case, we set S = 1/2 in the spin wave calculation. For 
0.22 << J 2 < 0.4, the pRVB state is lower in energy 
than the large-S estimate. 

(ii) 'Nematic Valence Bond Solid': a dimer state which 
breaks lattice rotational symmetry has been proposed as 
the ground state for J 2 > 0.4[19]. This state can also be 
viewed as the staggered dimer state on the honeycomb 
lattice. We plot the energy of this state obtained using 
bond operator mean-field theory keeping only quadratic 
terms in triplet operators [l^. For J 2 < 0.4, the pRVB 



state is lower energy than this dimer state. 

(iii) Variational Monte Carlo (VMC): recently, a VMC 
scheme using entangled-plaquette wavefunctions has esti- 
mated the ground state energy for the J\ — J 2 model [Tij. 
Around J 2 ~ 0.3, our energy estimate for the pRVB state 
compares well with the variational result. In fact, the 
pRVB energy is lower than the VMC result at J 2 = 0.3. 
Of course, our pRVB energy estimate may change when 
corrections such as triplet-triplet interactions, coupling 
to quintet states, etc., are included. While the VMC 
study does not find pRVB order, energy comparison in- 
dicates that pRVB is a plausible ground state. 

Comparing the energy of the pRVB state with the 
above candidates, we conclude that the ground state may 
have pRVB order for 0.25 < J 2 < 0.4. Next, we plot the 
spin gap as a function of J 2 / J\ - this is the energy of the 
lowest triplon state. We find that this lowest triplon state 
always occurs at the M points of the plaquette-Brillouin 
zone (see FigJHJ) . The spin gap reaches a maximum value 
of ~ 0.7Ji at J 2 ~ 0.225. Surprisingly, the spin gap does 
not close as J 2 is varied. We discuss its implications in 
the next section. 



V. DISCUSSION 

We have developed an effective theory for the proposed 
pRVB state in the spin-1/2 J\ — J 2 model on the hon- 
eycomb lattice. Using plaquette-operator mean-field the- 
ory, we have calculated the properties of this state and 
the spectrum of low- lying excitations above it. 

An important outcome of our calculation is to clarify 
the nature of the pRVB state. The low energy state is 
composed of a \[2> x y3 ordering of plaquettes in the '-' 
state. We begin our calculation from the limit of de- 
coupled hexagons, wherein the '-' state is clearly lower 
in energy for J 2 < 0.5. After including inter-hexagon 
coupling, we get a consistent theory of the pRVB state. 
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FIG. 8: The plaquette Brillouin Zone (BZ): The dashed line 
indicates the BZ of the underlying honeycomb lattice. The 
pRVB state occurs on a triangular lattice of plaquettes - the 
corresponding BZ is shown by solid lines using the extended 
zone scheme. The M points of the plaquette-triangular BZ are 
highlighted by blue/red circles - these are the wavevectors at 
which the lowest triplet excitations occur. 



plied magnetic field will lower the energy of triplet modes. 
At some critical field, the lowest triplet mode(s) will be- 
come gapless and undergo Bose condensation giving rise 
to magnetic order. In our case, the lowest triplet modes 
are triply degenerate as they occur at the three M points 
of the plaquette-BZ. Beyond the critical field, each mode 
can condense with an independent amplitude and phase 
leading to six degrees of freedom. This makes it difficult 
to predict the nature of the magnetic ordering beyond 
the critical field. At the level of these arguments, we can 
rule out field-induced Neel ordering as the lowest triplets 
do not occur at the T point. This picture may change 
once we account for defects in pRVB order as discussed 
below. 



B. Phase transitions 



Thus, our plaquette wavefunction changes sign upon ro- 
tation by 60°. Under rotation by angle <fi, the plaquette 
wavefunction transforms as e 1 ^ - having f-wave symme- 
try. 

To rule out s-wave ordering with condensation of the 
'+' state, we redid the plaquette operator calculation for 
this case. In the limit of disconnected hexagons, the '+' 
state is an excited state. Naively speaking, s-wave pRVB 
order is still possible if there is sufficient lowering of ki- 
netic energy. We find that this is not the case - the pla- 
quette operator description of s-wave pRVB order leads 
to unphysical results. We can thus rule out strong s- 
wave ordering - interaction corrections may modify this 
picture, but they should be unimportant when s-wave 
ordering is strong. 

The earliest study of plaquette order on the honeycomb 
lattice - using quantum dimer models - considered pla- 
quettes which are symmetric combinations [2(j- Howew- 
ever, for the J\ — J 2 model at hand, Refd indicates that 
the antisymmetric combination is favoured - the authors 
refer to this order as d-wave, while it is more precisely 
called f-wave. 



A. Excitations of pRVB state 

This paper presents the first systematic study of the 
pRVB state and its excitations. As shown in FigJTl the 
spin gap does not close within our theory. And the low- 
est lying triplet modes occur at the M points of the 
plaquette-BZ. In terms of the underlying honeycomb lat- 
tice, the lowest triplet state will occur at wavevectors 
shown in figUl Assuming that pRVB order occurs in the 
J\ — J2 model, numerical methods such as ED or DMRG 
can see this by calculating the spin-spin structure factor 
in the lowest triplet state. The structure factor should 
show peaks at the wavevectors shown in Figj8] Indeed, 
this can serve as a test for pRVB order. 

As is familiar from the study of dimer states (2lj. an ap- 



For the case of dimer-Neel transitions, it is well known 
that bond operator theory is remarkably sucessful in de- 
scribing the nature of the phase transition. Starting from 
the dimer state, bond operator theory provides a con- 
sistent theory of spin-1 triplon excitations which obey 
bosonic statistics. The transition to Neel order occurs via 
Bose-condensation of triplons when the spin gap closes. 

In our honeycomb lattice model, our plaquette oper- 
ator theory gives a consistent description of the pRVB 
state which is most stable around J2 ~ 0.25. ED calcula- 
tions suggest that a pRVB state is sandwiched between a 
Neel state(for J2 < 0.2) and a magnetic phase with stripe 
order (for J2 > 0.4). While the phase diagram definitely 
contains Neel order, the existence of a stripe magnetic 
phase is still open to question. It has been suggested 
that the ground state for J2 > 0.4 has dimer order which 
breaks lattice rotational symmetry [Til [l9j . 

Starting from the pRVB state, as J2 is decreased, we 
expect to see a Neel transition near J2 ~ 0.2. Naively, 
we may expect such a transition to be driven by conden- 
sation of triplons giving rise to a continuous Quantum 
Phase Transition (QPT) in the same universality class 
as the BEC transition in 2+1 dimensions. However, the 
spin gap obtained from plaquette operator theory does 
not close as J2 is lowered, thus ruling out a BEC transi- 
tion. This leaves us with two possibilities concerning the 
pRVB-Neel transition: (i) a first order phase transtion, 
or (ii) a deconfined QPT mediated by topological defects 
in pRVB order. 

In terms of conventional Landau theory, we generically 
expect a first order QPT as the Neel state and pRVB 
state break different symmetries. This is always a pos- 
sibility and can only be checked using precise numerical 
techniques. The second possibility, a deconfined QPT, is 
a very exciting prospect. So far, such a transition has not 
been observed in simple models involving Heisenberg-like 
couplings. 

A deconfined QPT would be driven by proliferation of 
defects in the pRVB order parameter - in this case, Z3 
vortices. As pointed out in Ref. 22], the core of such a 



8 



i + b I 



? 



\ * — IT 






1 i + a 1 


la — fel 



FIG. 9: A plaquette i and its nearest-neighbours, a and b are 
the primitive vectors of the \/3 x -^3 lattice of plaquettes. 

vortex will necessarily bind a dangling spin or a 'spinon'. 
The critical theory of a continuous pRVB-Neel transi- 
tion cannot be cast in terms of Landau-Ginzburg order 
parameters, but in terms of these spinful defects. We 
expect the theory to be very similar to the one outlined 
in Ref423] - at the critical point, the vortices become 
gapless and the Z 3 anisotropy may become irrelevant. 

Starting from the pRVB state, when we increase J2, 
there must be a transition into a state which breaks 
lattice rotational symmetry - a nematic Valence Bond 
Solid [HI or a magnetic phase with stripe order [l3j|. The 
spin gap does not close in our plaquette operator theory. 
As suggested by ED results |l3j, we surmise that this is 
a first order transition. Future numerical studies will be 
able to test this notion. 
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Appendix A: Matrix elements 

Fig. |9] shows a plaquette i and its six nearest neigh- 
bours. Let us first consider the nearest neighbour pair 
(i,i — b). This pair is identical to the that shown in 
Fig|4] if we identify (i) as A and (i = b) as B. Inter- 
plaquette coupling leads to hopping-like and pairing-like 
terms between triplon operators on the two sites. Below, 
we evaluate the coefficients h l ^~ h and d l £~ h . 

The generic inter-plaquette interaction process is indi- 
cated in Eq. [8] The hopping-like term arises the following 
matrix element. 

Ki-^iiiUi-bUH^iiiUi-b),}. (Ai) 

This matrix element can be evaluated numerically using 
two pieces of information: (i) the explicit eigenstates of 



the single hexagon problem and (ii) the explicit form of 
the inter-plaquette coupling given in Eq. [H In plaquette 
operator notation, this matrix element leads to the pro- 
cess b] ,b\ . bi _i>. ? ... As the '-' bosons are condensed 
with amplitude p, we rewrite this as 

H hoppmg ~ hit'VblA-h^ (A2) 
Similarly, the coefficient of the pairing term is given by 

dtit 1 = - luH^iaui - S)_). (as) 

This gives rise to b) ,, b ' t bi _ b . t , which can be 
rewritten as 

H pairmg ~ d%- l p 2 bl u bl_ iy (A4) 

Having evaluated these terms for one nearest neigh- 
bour pair, we can evaluate the matrix elements for all 
other pairs using the rotation operator. Let us consider 
three nearest neighbour pairs in Fig. |HJ (i,i — b), (i,i + a) 
and + a + b). Clearly, these three pairs are related 
by R, rotation by 60°, about the center of the hexagon 
i. Formally, in terms of the inter-plaquette couplings, we 
may write 

H^+a = //, , ,/,'. (A5) 

The rotation operator R maps the plaquette (i — b) to 
(i + a). In addition, it rotates the hexagon 60° about the 
center. 

Using this relation, we can obtain the hopping matrix 
element h^l +a : 

h^ +& = ((i) u (i + a)_\H^ l+ a\(i)-(t + a) v ) 

= ((»)„(» + a)_|#jj. ( ._ s A|(i)_(i + a) v ) 

= r u r:m a (i-bUH i ._ 6 \{i)_(i-b) v ) 

= r^;^, (A6) 

where r u and r v are the rotation-eigenvalues of the single- 
hexagon states ipu and ip v respectively. 

The coefficient of the pairing-like term can be found 
likewise: 

<S +a = ((*)«(< + a) v \H iti+& \(i)_(i + a)_> 
= ((»)«(* + a) v \R)H^ b R\(i)_{i + a)_) 
= r u r v (r*_) 2 ((i) u (i - b) v \H.^ t \(%)-(i - 

= r u r v {r*_fd^-\ (A7) 

Similarly, by applying one more rotation, we can ob- 
tain 

ii.i+a+b r 2 ( r *) 2 h i ' i ~' 1 _ „2 2 / * \4ji,i-S 

IL UV ' U\' V ) 11 uv j u ud 'u'vv—J u uv 
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